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The response of the system, consisting of two kinds of opposite-charged fermions and their bound states (hydrogen- like 
atoms), to the perturbation by the external electromagnetic field in low particle kinetic energies region is studied. 
Expressions for Green functions that describe the system response to the external electromagnetic field and take into 
account the presence of particle bound states (atoms) are found. Macroscopic parameters of the system, such as 
conductivity, permittivity and magnetic permeability in terms of these Green functions are introduced. As an example, 
the perturbation of the ideal hydrogen-like plasma by the external electromagnetic field in low temperature region is 
considered. Such approach also enables to study the propagation properties of the signal, tuned up to the transition 
between two hyperfine ground state levels of alkali atoms that are considered in Bose-Einstein condensation (BEC) 
state. It is shown that the signal can propagate in such system with rather small energy loss. Such fact allows to 
introduce the group velocity concept and to study the slowing down conditions for the microwave signal that propagates 
in BEC. 
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1 Introduction 

In the process of describing a behavior of many-particle 
systems a class of problems appear, that concerned 
with the system response to the perturbing action 
of the external, in particular, electromagnetic field. 
Widespread approach to solving such kind of problems 
is based on using the Green functions formalism (see 
in that case e.g. Q]). 

As well known, the most convenient method 
for physical processes describing in quantum many- 
particle theory is the second quantization method. Be- 
cause of that in framework of the second quantization 
it is the most simple to formulate an approach for de- 
scribing the system response to the perturbation by the 
external field, that is based on using Green functions. 
But if we try to realize such an approach, we can meet 
an essential difficulty, connected with the possibility of 
the particle bound states existence. 

Really, the key role of the second quantization 
method consists in the introduction of creation and 
annihilation operators of particles in a certain quan- 
tum state. The operators of physical quantities are 
constructed in terms of creation and annihilation op- 
erators. Such a description of quantum many-particle 
systems implies the particles to be elementary (not 
consisting of other particles). Moreover, it is abso- 
lutely accurate despite of the possible existence of 
compound particles. Since the interactions between 



particles may lead to the formation of bound states, 
the standard second quantization method becomes too 
cumbersome. For this reason the construction of an 
approximate quantum-mechanical theory for many- 
particle systems consisting of elementary particles and 
their bound states represents an actual problem. In 
this theory it is necessary to introduce the creation and 
annihilation operators of bound states as the operators 
of elementary objects (not compound). In addition it 
must preserve the required information concerning in- 
ternal degrees of freedom for bound states. 

Such an approach is realized in [2] . In this work the 
possibility of constructing such theory is demonstrated 
for a system, that consists of two kinds of fermions, 
assuming that bound states (atoms or molecules) are 
formed by particles of two different kinds. The choice 
of such model is not dictated by the principle diffi- 
culties but by the desire to simplify calculations and 
obtain the visual results. In framework of this model a 
method of constructing the creation and annihilation 
operators of the bound state as a compound object is 
given. The substantiation of the conversion from the 
description of atoms as compound objects to elemen- 
tary objects with the ordinary creation and annihila- 
tion Bose-operators is given. Such substantiation is 
considered in low-energy approximation in which the 
binding energy of a compound particle is much greater 
than its kinetic energy. In terms of the creation and 
annihilation operators of fermions and bosons (as ele- 



mentary objects) a scheme for physical quantities oper- 
ators construction is formulated. Explicit expressions 
for the operators of principal physical quantities, such 
as density and charge density momentum and current 
density system Hamiltonian, are found. The Maxwell- 
Lorentz system of equations, describing the interaction 
between electromagnetic field and matter, that may 
consist also of neutral "atoms" (low-energy quantum 
electrodynamic equations), is found. 

In the present work we use this system of equations 
to study in the framework of Green functions formal- 
ism the response of the system with bound states to 
the perturbation by the external electromagnetic field. 
An essentially new moment in these considerations is 
the next circumstance. When we describe the system 
response to the perturbing action of the external elec- 
tromagnetic field the approximate formulation of the 
second quantization method proposed in [5] gives us 
the possibility to account for the neutral bound states 
in sufficiently simple way. 

2 Quantum electrodynamic 
equations for the low- 
temperature hydrogen-like 
plasma 

The studied quantum-electrodynamic system, consist- 
ing of fermions of two different kinds and their bound 
states, in low-energy region, in fact, can be considered 
as a low-temperature hydrogen-like plasma. Before we 
turn to the description of such plasma response to the 
external electromagnetic field, let us obtain the main 
equations that describe an evolution of such system. 
Taking into account the interaction between radiation 
and matter the system's Hamiltonian H(t), according 
to [2], can be written as 

H(t)=Ho+H int + V(t), H =H f +H p , (1) 

where 



(2) 



k.A 



is the Hamiltonian for free photons (u>k is the frequency 
of photon with wave number k, C^Ai are the cre- 
ation and annihilation operators of photon with wave 
number k and polarization A). 

The value H p in the equation (flj is the Hamiltonian 
for free particles (free fermions and their bound states) 
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where x^( x ), Xj( x ) (j = 1,2) are the creation and 
annihilation operators of a free fermion of j kind and 
mass m,j at the point x; ?)+ (X) , f] a (X) are the creation 
and annihilation operators of bound states of two dif- 
ferent fermions (" hydrogen- like atoms" ) with the quan- 
tum numbers a at the point X; e a is the energy of an 
atom at the level with the quantum numbers a. 

The Hamiltonian TLi n t in the equation (fT]) describes 
the interaction between particles 



n 



int 



n_j 1 _|_ a_/2 



o_/3 



(4) 



where 



= J dx 1 dx 2 dy(p + (x 2 ,y)0(x 2 ,y) 

x {(i/u(xi - x 2 ) + i*2i (x x - y))Xi~( x i)xi(xi) ( 5 ) 
+ (1^22 (xi - y) + via (xi - x 2 ))x 2 f (xi)x2(xi)}, 



dx 1 dx 2 dy 1 dy 2 + (x 1 ,y 1 ) 

x ^ + (x 2 ,y 2 )^(x2,y 2 )^(xi,y 1 ) (6) 
x {fu(xi - x 2 ) + ^ 22 (y x - y 2 ) 
+ i/ia (xi - y 2 ) + ^21 (Yi - x 2 )}, 



= g / dx x dx 2 
x {i/u(xi - x 2 )Xi (xi)Xi f (x2)Xi(x2)xi( x i) (7) 

+ 1/22 (Xl - X 2 )x2"(Xl)X2 (X2)X2(X 2 )X2(X1) 

+ 2i/i 2 (xi - x 2 )xi" (xi)x2 (x 2 )x2(x 2 )xi(xi)}. 

In equations ([5]) — (0) the operators <^ + (xi,x 2 ), 
<p(xi,x 2 ) are related to the creation 77+ (X) and an- 
nihilation operators f) a (X) of atoms in quantum state 
a by expressions 

^+( Xl ,x 2 ) = ]r<(x)7)+(x), 

<p{xi,x 2 ) = ^2(p a (x)fj a (X), (g) 

a 

mixi + m 2 x 2 



x = xi — x 2 , X 
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where ip a (x) is the wave function of the bound state 
and fij(x — y), i,j = 1,2 is the potential energy of 
Coulomb interaction 



i^(x-y) 



(9) 



(ej is the electric charge of a fermion of kind i). 

In that way, the Hamiltonian Hj ni corresponds to 
scattering of particles of the first and second kinds 
by bound states, the Hamiltonian nf nt corresponds to 
scattering of bound states by bound states, the Hamil- 
tonian li? int corresponds to scattering of particles of the 
first and second kinds by particles of the same kinds. 



And, finally, the operator V(t) in (JTJ) represents the 
Hamiltonian that describes the interaction of particles 
with the electromagnetic field 



1 



V(t) =-- I dxA(x,i)J(x,i) 
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^|dxA 2 (x,;)]T-^(x) (10) 

J dxip ( - e \x, i)o-(x), o-(x) = y^o-j(x). 



In this expression we have taken into account an inter- 
action of particles with the external electromagnetic 
field A (e) (x,i), ^ (e) (x,i) (ip( e )(x,i) is the scalar po- 
tential of the external electromagnetic field) and the 
quantum electromagnetic field, that is described by the 
potential a(x) (Coulomb's gauge): 



A(x,t) = a(x)+A (e) (x,i), 



(11) 



where A^ (x, t) is the vector potential of the external 
electromagnetic field and a(x) is the quantum electro- 
magnetic field operator, that is defined by the expres- 
sion 
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(V is the system volume, ekA is the photon polarization 
vector). 

The charge density operators <Xj (x) for particles of i 
kind (see (jXOj) ) are connected with the density opera- 
tors /3j(x) (see [2]) 

o-j(x) = e i p i (x), 

Mx) = tf(x)*i(x) + f d yf dY 

x<5(x-Y-^y)^+(y,Y)^(y,Y), (12) 
^x) - xj(x)x 2 (x) + J dy J dY 

x5(x-Y-^y)£+(y,Y)£(y,Y), 

and, as easy to see, in the equation (fT2")l we have also 
taken into account a contribution that had been made 
by charged particles, that are represented in the bound 
states (see ®). The current density operator J(x, t) 
in the formula (JTOj) can be also expressed in terms of 
the creation and annihilation operators 



J(x,i) = -A(x,t)Y]-^a i (x)+j(x) ) 

z — ' TO, 
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where the momentum density operators 7Tj(x) are de- 
fined by expressions 
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Using formulas (fl"2"l) ~ (fl"4]) we can write the expres- 
sions for the current and charge density operators in 
more suitable way: 



o-(x) = ^£ a (x), 

a 

j(x)=^j„(x), a = 0,1, 2, 
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where 

CTj(x) = e,x+(x)xi(x), i = 1,2, 
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o-o(x) 
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As is easy to see, the operators <xo(x), j (x) in these 
expressions define the bound states contribution to the 
charge and current densities. 

In the momentum representation 



^(x) = -L^e^fy Q (p) 



expressions (jT6|) according to ([8]) will have the next 
form: 
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+ " 2M 
where (see ([5])) 
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It is significant to note, that the Hamiltonian for free 
particles H. p (see ([lj, ([3])) in the momentum represen- 
tation can be written as 



Hp — H\ p + H 2p + Hq p , 



= z2 £i(p)at p a ip , i = l,2, 
p 

Hop ^EE^fp^fp)^"^' 



(19) 



a p 



£i(p) = P 2 /2nii, e a (p) = e a + p 2 /2M, 

where e a is the energy of the atomic level with quantum 
numbers a, M is the bound state mass, M = (mi + 
m 2 ). 

The Maxwell equations for our system according to 



[2] can be written in the following form 

— = -crotE, divE = 4tt(<7 + cr (e) ), 
ot 

— = crotH-47r(J + J (e) ), divH = 0, 
at 



(20) 



where operators <r, J are still defined by the expressions 
CEO]), Cll), CE3]) and values cr (e) , J (e) are the external 
current and charge densities. The electric E and mag- 
netic H field intensity operators in terms of the scalar 
and vector potentials can be expressed as (see pQ, [2] 
and also equations (fTU|) . (fTTjl ) 



H = rot A, 



E 



1 dA d 

c dt dx 
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dx 



, a(x') \ (21) 



Note that in deriving the electrodynamic equations we 
used the Coulomb's gauge. 

3 The system response to the 
perturbation by the external 
electromagnetic field 

In this section to study the system response to the 
perturbing action of the external electromagnetic field 
we will follow the principles that have been stated in 
[1]. Let us consider a system that at some moment 
of time t is characterized by statistical operator p(t). 
Considering that the Hamiltonian of interaction V(t) is 
linear in respect to the external field and assuming that 
it is small in comparison with the Hamiltonian H = 
Ho + Hint (see ([TJ), we can develop the perturbation 
theory over the week interaction. In accordance with 
[T] the mean value of an arbitrary quasilocal operator 
a(x) in linear approach for such system can be written 
as 



Spp(t)a(x) = Spwa(0) +a F (x, t), 



a F (x,t) 



dt' / dx' 



(22) 



xG^ix-x'^-mix'^'), 
where w is the Hibbs distribution operator 

w = exp{ft - (3{H - mNi - H2N2)}, 



(23) 



[3 = 1/T is the reciprocal temperature, N\, iV 2 are the 
density operators of all fermions of the first and second 
kind (including fermions in bound states, see (|12p ) 



dxpi(x), N 2 = / dxp 2 (x) (24) 



Ni 



and /ii, /i 2 are the chemical potentials of fermions of 
the first and second kind. The thermodynamic param- 
eters (3, fii, p 2 can be found from the relations 

SpwH = H, SpwTVi = N x , SpwN 2 = N 2 , (25) 



and the thermodynamic potential f2 dependence on 
thermodynamic parameters is defined by the expres- 
sion 

Spit; = 1. 

In the formula (f!2"2"|) Fi(x,t) are the quantities, that 
define the external field and £j (x) are quasilocal opera- 
tors, related to our system (see also [I]); the summation 
convention is assumed for the repeated index i. 

And, finally, the quantity G^(x — x',t — t') in the 
expression (|22l) is the two-time retarded Green function 
(note, that "tilde" over operators means that they are 
taken in the Heisenberg representation) 



(26) 



G { a f i (x-x',t-t') = -i9(t-t') 

xSpw[a(x,t),ii(x',t% 
where 9(t) is Heaviside function 

*(t) = ( M>0 ' 
w [0, t < 0. 

Going over to Fourier transforms of values a F , Fi 
a F (x,t) = J dkd W e-^- kx )a F (k,c), 

Fi(x,t) = J dkckve-^-^F^u;) 

one obtains 

a F (k,u) = G ( +l(k,oj)F i (k,u ) ), (27) 

where 

G%>Ql,u)= dt / dxe* w - kx )Gg(x,t). (28) 



It is significant, that in terms of the Fourier trans- 
forms of the introduced quantities we can express also 
the energy, transferred from field to matter. If to as- 
sume, that field is acting only for a limited period of 
time, the total energy Q, received by matter, is given 
by the expression pQ 



Q = 



(27T) 



dio / dkuj 



xFi(-k, -u)G$l (k, wjF^k, w). 



(29) 



Now we can apply these expressions for studying the 
response of the system, consisting of two kinds of op- 
positely charged fermions and their bound states. To 
make use of the Green function method, that was de- 
scribed above (see (|22|) -(f28 | ). it is more convenient to 
represent the system Hamiltonian, that is defined by 
formulas (fT])- (fT6|) . as 

H(t)=H + V^(t), H = H a +H mt + V, (30) 



where Ho and Hint are given by the formulas ([I])-®, 
V is defined by the expression (see also (fT0|)-(fT4l)) 



V = -- I dxa(x,i)j(x,i) 



If 2 
^ydxa 2 (x,t)£^(x), 

J i=i 1 



(31) 



and the Hamiltonian V^ e '(t) describes the system in- 
teraction with the external electromagnetic field 

V^(t) = ~ J dxA (e Hx,<)j(x) 

+ ^2 /dxA( e >(x,i) 2 ]T-^(x) (32) 



To get the Maxwell equations for the electromagnetic 
field in medium, it is necessary to average the equations 
(|20p with the system statistical operator containing the 
information both about medium and electromagnetic 
field. To this end we will define the mean values of 
the electromagnetic fields E(x, i), H(x, t), acting in the 
matter 

E(x,i) = Spp(t)E(x,t), H(x,i) = Spp(t)H(x,t), 

(33) 

and also the induced charge and current averages (see 

J(x, t) = Spp(t) J(x, t), cr(x, t) = Spp(t)<r(x). (34) 

The equations (|20|) . averaged in accordance with the 
formulas (|33 | - ([34|) . bring us to the Maxwell-Lorentz 
equations for the average fields in the matter 



<9H 

~dt 
<9E 



-crotE, divE = An{a + cr (e) ), 
crotH-47r(J + J (e) ), divH = 0. 



(35) 



where quantities <j^ e \ 3^ still represent the extrinsic 
charge and current densities. 

The next problem is to find the charge <r(x,i) and 
current J(x, t) densities induced by the external field. 
Calculating these quantities under assumption of week 
interaction between the system and the external field 
we will use the equations ([2"2")) - (f2"g]) . considering the 
potentials A^(x, t), <^ e )(x, t) as Fi(x,t), and <r(x) or 
J(x, t) as a quasilocal operator a(x). As a result one 



gets 



5(x, t) = ^2<r a + J dt' J d 3 x' 

— OO 

-Gi +) (x-x',t-OiA[ e) (x',<') 

+GW(x-x',t-t')^ e) (x',<')] , 

J fc (x,i) = -i4 e) (x,i)]T^a Q 
c ^ m 

a 

oo 

+ J dt'Jd 3 x' -Gj^x-^.i-O^V.O 



+Gj c +) (x-x',i-<> (e) (x',0 

(36) 

where cr a = Spw Q (0), a = 1,2,0 (see (fl7|)). w is the 
Hibbs statistical operator (|23|) . It is significant to note, 
that for quasincutral systems, where the number of 
fermions are equal (Ni = N2, see (j2"4")l ) and their abso- 
lute charge values are also equal (|ei| = 1 62 1 ) Y^cr a = 0. 

a 

The retarded charge and current Green functions, 
that are included in the expression (|36|) , are determined 
in accordance with the formula (|26|) (see also pQ): 

G (+) (x,i) = -i6(t)Spw[a(x,t),a(0)}, 

Gi +) (x,t)--j0(t)Sp«;[j fc (x,t), ( 7(O)], 

-(+) - ( 37 ) 

Gi +j (x,t) = -i0(*)Spu>[<7(x,t),j fc (O)], 

G<+ } (x,i) = -*fl(i)Sp«;^(x,t)J / (0)]. 

As the charge and current density operators of particles 
of different kinds (see (fTT|0 commute with each other 

[<7a,<7fc] = [<TaJ b ] = [j a ,j b ] =0, a ^ b, a, b = 1,2,0, 

then, according to the equations (fT5)l . (JTTJ), (|3T)l . the 
contribution of different kinds of particles to Green 
functions will be additive 

GW(x,i) = ^GW(x,t), 

a 

Gi +) (x,<) = -iO(t)Spw[a a (x,t),a a (0)], 
Gi +) (x,t) = ^Glt ) (x,i), 

a 

( X )<) = -^(0 S P w [^a(x,t), jafe(0)], 

Gi +) (x,t) = ^GW(x,i), 

a 

Gj+W) = -»0(i)Sp W [j Qfc (x,t), ( T Q (O)], 

a 

G$(x,i) = -<fi(t)Sp«;[i ofc (x,t),5- ol (0)]. 

With the help of direct calculations, following the 
method p], we can see, that also in the presence of 



(38) 



particle bound states the following correspondence be- 
tween Green functions (f3"T)) takes place 



Gi T '(x,i),Gi +) (x,t), 
dG< +) (x,<) 0G(+)(x,i) 



St 



0. 



SGW(x,t) , 0G< +) (x,t) 



(39) 



at 



E— ^(t)^-J(x)=0 

* ■ 771 ' ' ■' 



For the Green functions Fourier transforms (see (|28[)) 
these relations can be written as 

Gi +) (k,c) = Gi +) (k, w ), 
G| +) (k )W )fc i - W G(+»(k )tl) ) = ) (40) 
G 4 ( + } (k, u)kj - wG| +) (k,w) + £ -t« - 0. 



4 Green functions and macro- 
scopic characteristics of 
the ideal low-temperature 
hydrogen- like plasma 

If we neglect of all interactions between particles in the 
investigated system, it can be considered as an ideal 
hydrogen-like low-temperature plasma (we note, that 
kinetic energy of particles must be small in compari- 
son with the binding energy of compound particles). 
For an ideal hydrogen-like plasma the Green functions, 
that was introduced earlier, can be calculated exactly. 
To do it we must take into consideration, that neglect- 
ing the quantum fields presence, the Hamiltonian Ti. in 
the formula ([30]) must be interpreted as H p , see (pQ), 
((3]), (fl9|) . Taking into account this fact the Hcisenbcrg 
representation of charge and current density operators, 
that appear in the equations (|37p for the Green func- 
tions, is defined by expressions: 

a,(x,t) = ^]>>- 4X (P-P0 e -«(^(P)-^{p')) a + a4p , j 



P:P' 



Ji(x, t) = <7i(x, t), i = 1, 2, 

2m, i 



-p') p - rf (e»(p)-e/3(p')) 



P,P' a, 13 

x a a p(p - p')fi+(p)f) (p'), 
j ( X;t ) = ^^^ e -ix(p-p') e -«(«-(p)-«/>(p')) 



P,P' a, 13 
(p + p" 



2M CT ^P-P') 

+Iaf3{p-P')}fla (P)^(P'). 



(41) 



where quantities <7 Qi g(k), I Q ^(k) are given by the for- 
mulas (|T5| . If to substitute the operators (|4*Tj) in 
and to do some calculations, we will come to the fol- 
lowing expressions for the Fourier transforms of scalar 
Green functions (see (12810: 



r«fk , A - e ' V A(p-k)-/i( P ) 
tk,u, "y^ ei (p)- El ( P -k) + u + ! o' 



r (+) (v , a _ e l V- / 2 (p-k)-/ 2 (p) 
2 ( ' j " V ^£ 2 (p)-£ 2 (p-k)+c + ? :0' 

G< +) (k,c) = I^^^k^-k) 



P ",/3 



/q( P -k)-^(p) 

e Q (p) - e /3 (p - k) + uj + iO' 



(42) 



Similar expressions for the vector Green functions have 
the form: 



G\V(k,oj) = _^_y^(2p-k) 
11 y ' ' 2 mi V ^ V y ' 



G^(k^) 



/i(p-k)-/!( P ) 
£i (p) - £i(p — k) + oj + iO ' 



2m 2 V 



/ 2 (p-k)-/ 2 (p) 
£2 (p) - £ 2 (p - k) + w + iO ' 



G, 



(+) 

0/ 



(2 P - k) 

2M 



CT Q/3 (k) +I a/3 (k) 



P a, (3 

<W~ k ) [/q(p ^ k) - fp(p)] 

e a (p) - e / a(p - k) + w + iO ' 



(43) 



And, finally, the tensor Green functions for the inves- 
tigated system is given by expressions: 

Gffi (k, «) = £ (2P ~ k) ; (2 P - k) s 

1 p 

x /i(p-k)-/i( P ) 
ffi(p) - ei(p — k) + w + i0' 



4m2V^ 



5^(2p-k),(2p-k) 4 



/ 2 (p-k)-/ 2 (p) 
£ 2 (p) - e 2 (p - k) +cj + iO' 



P a, p 



(2 P - k) 

2M 



o-a^(k) +I a/3 (k) 



(2 P - k) 

2M 



CT/3a(-k) + I / 9„(— k) 



/a(p-k)-/ff( P ) 

£q(p) - £/3(P - k) + w + iO' 



(44) 



In the expressions (fl2"|) - (|4"4")) we have introduced the 
distribution functions for free fermions of the first fi (p) 
and second / 2 (p) kind, and also the distribution func- 
tions / Q (p) for hydrogen- like atoms (bound states) 
with the set of quantum numbers a 



(45) 



(46) 



/ I (p) = {ex P [(e l (p)-^)/r] + l}- 1 , 
/ a (p) = {exp[( £Q (p) - /z Q )/T] - I}" 1 

in accordance with the equations 

Spwaf p a ip , = 6 PtP >fi(p), i = 1,2, 
Spwf?+(p)?) / 3(p / ) = 6 at p6 PtP >f a (p). 



The particle energies £1,2 (p), £ Q (p) in formulas f4*2|l 
(|45ll are given by the expressions (fl9|) . and values 5 aj( g, 
(5p, p ' in (|46|) are the Kronecker symbols. 

The particular feature of the obtained Green func- 
tions is that the contribution of bound states in the pro- 
cesses under consideration now is taken into account. 

The expressions for Green functions that have been 
found allow us to get an expressions for the matter 
macroscopic parameters, such as conductivity, permit- 
tivity and magnetic permeability. To this end we will 
also use the method described in [1] . 

In accordance with the formulas (f35|) . |36|) the next 
relation between Fourier transforms takes place 



Ji(k,w) =a l {k,oj)k. 



kE (e) (k,w) 



1 



k 2 

[k,E< e )(k,u;)],k 



(47) 



SF(k,w) =-5'(k,Lj)kE (e) (k,w), 



where 



^(k.w) 
^*(k,w) 



'k 2 

i 



G (+) (k,w), 



a 



(48) 



2 ,4 +) ( k - w ) 



It is clear from equations (|4T|) that the quantities a J 
and a , expressed in terms of Green functions according 
to formula (f48|) , define the longitudinal and transversal 
current density components. They are usually inter- 
preted as outer conductivity coefficients in contrast to 
inner (or true) longitudinal a 1 or transversal cr* con- 
ductivity coefficients, that will be defined below. Note 
that according to (|38|) these coefficients are also addi- 
tive quantities 



^(k,c) = E^'(k, 



In terms of the introduced outer conductivity coef- 
ficients (7 and a 1 we can express the energy, absorbed 



from the external field sources (see 



Q.k = jtt 4lm-E* (e) (k,c)G+(k^)E( e )(k, W ). 
[Zirp lu J 

From this expression, in accordance with the formulas 
(ETb-dHh one obtains 



(2tt)4 



Re{a / (k, W )|E[ | e) (k, W )| 2 
+* t (k,u J )\E ( l\k,u J )\ 2 }. 



(49) 



In terms of these quantities (a and tr ) the expres- 
sions for permittivity and magnetic permeability can 
be also defined. The relation between the permittivity 
and outer conductivity (see [1] ) is given by the formula: 



e= 1 



47TO-' 



ILO 



(50) 



From this, according to the expression 
e-^w) = l + -jG«(k,< 

It is more convenient to express the magnetic perme- 
ability in terms of inner conductivity coefficients a 1 and 
a* ' 

/i- 1 (k^) = l + ig^( ( 7 Z - ( 7*), (51) 

that are connected with the outer conductivity coeffi- 
cients a 1 and a 1 by the relations: 

i -i t &t 

a = ecr, a = 



fc 2 r 2x _l ' 



47TCT* , 

1+-^- I 1 

ILO 

In accordance with the expression (|49p. their rela- 
tions with the Green functions (l42l), (l44l) take the form: 



a\k,io) = 



fc 2 +47rG(+)(k,w)' 



A(k,w) 



CT ^ k,W ' ) iw (w 2 -/c 2 c 2 )+47rA(k,w)' 



(52) 

So, we have defined the main macroscopic charac- 
teristics of the ideal hydrogen-like plasma in low tem- 
perature region. These characteristics allow us to solve 
a number of applied problems for our system. Let us 
demonstrate it on few examples. 

Using the developed theory it is not difficult to find 
the permittivity of an ideal gas of hydrogen- like (alkali) 
atoms at low temperatures. According to the expres- 
sions (j42|) . (|50|) in neglect of free fermions contribution 
one gets 



e \k,Lu) = 1 + Tj^E^W 17 ^^^ 



k 2 V 

P a,P 

, /q( P -k)-^(p) 

N e Q (p) -E /3 (p-k) +w + i0' 



(53) 



As good known, at extremely low temperatures the 
Bose- Einstein condensate (BEC) of alkali atoms can be 
formed. At the temperatures much lower the critical 
point temperature (T <C T , see e.g. [I]) the bound 
states distribution functions f a (p) are proportional to 
the Dirac delta-function <5(p). Therefore, according to 
the expressions p^|) . (|53[) . after integration over mo- 
mentum p the expression for permittivity of the stud- 
ied gas in BEC state (T <C Tq) take the form: 

e _1 (k,u;) « 1 + 7ip^^W^(-lt) 



a,P 



lu + Ae a/3 - £fe + i7 Q /3 

^ 

to + Ae a/3 + e k + i~i afL 



(54) 



where v a is the density of condensed atoms in the quan- 
tum state a, Sk = k 2 /2M and quantities a a pik) are 
still defined by the formula (fl~8|) . Note that due to the 
damping processes in real systems we also introduced 
linewidth ~/ a p, concerned with the transition probabil- 
ity from the state a to the state (3. 

As it easy to see, in the expression (l54l) at frequencies 
that are close to the energy interval Ae Q ^ (Ae a/ 3 = e Q — 
sp) some peculiarities appear. In fact, such behavior 
can strongly reflect on the dispersion characteristics of 
the studied gas. It is a very interesting question, and 
we shall return to it in the section [5] 

Basing on the developed theory we can also find the 
energy, that is dissipated by a charged particle when it 
passes through hydrogen-like plasma at low tempera- 
ture (see in that case, e.g. PQ). In the case of a small 
dissipation the particle movement can be considered as 
uniform. Thus, the particle current density (the exter- 
nal current density in medium, see (|20p ) will be defined 
by the formula 

J (e) (x,f) = zev<5(x- vt), 

where ze is the particle charge and v is the particle 
velocity. It is easy to see, that the Fourier transform 
of the current density is given by the expression 



J (e) (k,w) = 27rzev(5(w-kv). 



(55) 



Next we will use the expression (|4"9"]l for the energy, ab- 
sorbed in the matter from external field sources. In this 
expression, using (|47|) . the Fourier transforms of the 
longitudinal and transversal components of the exter- 
nal field can be expressed in terms of longitudinal and 
transversal components of the current particle density 
(|53|) . If we do the necessary calculations, we will come 
to the following expression for the energy dS^, that 
was dissipated by the charged particle per unit time 
in the frequency dio and the wave vector dk intervals, 
when it passes through the hydrogen-like plasma: 



<7k^ = 



T 



(g) *(o,-kv) W 

fc2 x-i (56) 



\ c 2 e/i/ \ c 2 /i 

where T is the particle time of flight. To get the ex- 
pression (|56p it is necessary to use the formula 

<5 2 (w-kv) = — 5(co-kv). 

The total dissipated particle energy £ per unit length 
can be found by integrating the expression (|56p over lu 
and k 

d£ 1 f , 

— = -- / dcod 6 kq kLJ . (57) 

It is easy to see, that the main contribution in this 
integral comes from poles of the integrand (see (|56[> ) 

lu 2 

e(k, cj) = 0, -^-e(k, wWk, w) - fc 2 = 0. (58) 
c 2 

The formulas (|56]) - ([55)) are similar to the expressions, 
that are given in pQ, however, they take an account of 
particle bound states (atoms) to all processes that take 
place in our system (see p|-([IIJ), ([4"2 j) -(|4"6 ]l ). Note 
also that the expressions (|58|) represent the dispersion 
relations for free waves, that can spread in the studied 
system. 

5 Light delay phenomenon for 
the two- level system in BEC 
state 

In the previous section the expression for the permit- 
tivity of the system in BEC state (see Eq. ([54]) ) was 
found. For the system with the frequency of the ex- 
ternal laser field tuned up to the difference between 
two defined levels (marked below by indexes 1 and 2, 
to rs Ae 2 i) it can be written in more suitable form: 



e-\k,u))^l 



fflg2|g!2(k)| 2 [Vl - V 2 ) 
2lT 2 k 2 



(59) 



5ui + «7 

Here gi is the degeneracy order of i level, Slu — lu— Ae 2 i 
is the laser detuning, 7 = 7x2 is linewidth related to 
the transition probability from the upper to the lower 
state. Note that we also neglected the term £/., the 
substantiation of such operation will be discussed be- 
low. 

To study the propagation properties it is more con- 
venient to turn to the refractive index and damping 
factor quantities. To do it we set the magnetic perme- 
ability that enters the dispersion relation (|58|) close to 
unity fi(k, to) = 1. In that case the refractive index and 
damping factor are concerned with the permittivity by 
well-known relations: 



y/e' + Ve' 2 + e" 2 
V2 



7T~ 



C II2 



(60) 



Here e' and e" are the real and imaginary part of the 
permittivity: 

e(k,w)=e / (k,w) + ie"(k,a;), 

which can be found from the expression (|59p by taking 
the real and imaginary part: 



5u)(5w + a) + 7 2 
JEu + a) 2 + 7 2 ' 



(5u> + a) 2 + 7 ; 



r, (61) 



where 



a(k) = (vi - v 2 ) 



gig2ki2(k)| ; 
2vr 2 fc 2 



(62) 



Now one can find the dependence of the group ve- 
locity on the system parameters: 



(63) 



n + Lu(dn/duj) 



As it easy to see, the ultraslow light phenomenon can 
be observed in case when matter has a rather strong 
dispersion, or, according to the formula 



uj(dn/du)) > 1. 



(64) 



It is known that to use the group velocity concept the 
energy dissipation must be rather small, or, in other 
words, the next condition must take place: 



k'l « |e A 



(65) 



If e! ~ 1, then according to the expression (|60| and re- 
lation (f6"5)) (dn/dui) « 0.5(de' /duo). The partial deriva- 
tive of the permittivity real part can be found from the 
formula (IBTl): 



de' _ a [{Slu + a) 2 - 7 2 



du; [(Slu + a) 2 + 7 2 ] 



(66) 



Using Eqs. ([63 j) - ([66|) for the slowed signal group veloc- 
ity finally we get: 



((^ + a) 2 +7 2 ) : 



Vg ra 2C au [(Slu + a) 2 - 7 2 ] 



(67) 



Now let us consider the cases in which the light delay 
phenomenon for the two-level system in BEC state can 
be observed. To do it we shall go to the limit Slu — > 
0. In that case according to Eq. (|6T]) the real and 
imaginary parts of permittivity will have the following 
limits: 



lim e = 



T 



7 2 + a 2 ' 



lim e" 



1 a 



7 2 + a 2 



Now one can see that the condition for the dissipation 
smallness (I55|) is equivalent to the following condition: 



a 

— < 1, 

7 



Note, that it is necessary also to add the slowing down 
(strong dispersion) condition (|64[) . which in this case 
can be written as: 



c 
v„ 



Ae 2 i 1 ^ > 1, 



As for the defined system with fixed energy spectrum 
structure the one parameter that can be changed is oc- 
cupation difference Sv = (y\ — 1/2), which is included 
in the parameter a (see dM}), thus, basing on these 
relations we get the expression that characterizes the 
region where the mentioned phenomenon can be ob- 
served: 

7 \a\ 



< — < 1, 
Ae 2 i 7 



(68) 



Let us demonstrate that such region can exist on the 
example of cesium hyperfine structure ground state lev- 
els. The choice of hyperfine levels is stimulated by their 
stability and pumping capability (see more in that case 
in Ref. [3]). Note that for such levels the dipole tran- 
sitions are forbidden, thus the transitions come from 
the higher order effects that reflects on the extremely 
little values of linewidths. It will be shown below that 
such fact gives the opportunity for signal to propagate 
with small loss of energy. 

It should be mentioned that the description analog- 
ically can be spread to the other hydrogen-like atoms 
and other type of levels. For example, it can be used 
for description the experiments with ultraslow light in 
BEC of sodium atoms [3], where for the pulse slowing 
the three-level system had been used. 

As good known, alkali metals ( 133 Cs in particular) 
in the ground state do not have the dipole moment, 
thus the charge density matrix element oyi (see defini- 
tion (fT8")Q must be expanded to the second order over 
(ky) <1. As a result, one gets: 



0-12 (fc) ~ g( fcr, C)) 2 , 



(69) 



where r$ is the atomic radius (for the cesium ground 



state ;•() 



2,6 x 1CT 8 cm M 



is the electron 



charge. Taking g\ = 7, g2 — 9, the lincwidth 7 ss 
3, 8 x 10~ 21 cV (or 10~ 6 Hz in frequency units, as in 
"cold atomic clocks" experiments [6]), k ~ (Ae 2 i/c), 
where A^2i ~ 3, 8 x 10 -5 eV (microwaves with fre- 
quency 9,2 GHz) and basing on the expressions (f62|) . 
(|68p one can find the next region for the occupation dif- 
ference when light delay phenomenon can be observed: 



10" 



< \vi -v 2 \ < 3 X 10 13 cm 



(70) 



In the experiments for the BEC regime reaching the 
cesium atoms with peak density of v = 7 x 10 10 cm~ 3 
was kept in the trap [7], so, according to the inequality 
(|70| . one can conclude that also for such experiments 
the microwave pulse slowing down phenomenon can be 
observed. Note that the effect becomes greater with the 
density difference increasing until it reaches the upper 



limit of the expression (I70j) . when damping effects in 
the system are great. 

Also we must note the following. As it easy to see 
from Eq. (j6"T|) in the limit 5uj — > 0, the sign of the 
group velocity v g depends directly on the sign of the 
quantity a, that in turn depends on the sign of the 
difference (v\ — v-i). In other words, it depends whether 
population is normal or inverse. 

In case of normal population (y\ > v-i) the group 
velocity of signal is negative. It is traditionally consid- 
ered that the group velocity for the transparent matter 
is positive. But here one can conclude that due to the 
relation |68|) . the signal can propagate in the system 
with rather small dissipation (i.e., in fact, the matter 
is transparent) and rather slow velocity. Let us note 
that an observing of electromagnetic pulses with neg- 
ative group velocity is not so abnormal. The existence 
of such kind of phenomena for physical systems when 
the wave frequency is close to atomic (or molecular) 
resonances was pointed out and studied in many works 
(both theoretical, see e.g. [8], and experimental [91IT0]). 
In case of inverse population (y\ < 1/2) more "normal" 
situation takes place because the group velocity of the 
slowed pulse is positive. 

We stress that such rather unusual phenomenon oc- 
curs due to the unique property of hyperfine splittcd 
ground state levels (|69|) . One can show that for the 
two-level system with allowed dipole transitions the 
signal will not propagate due to large absorbtion. The 
most obvious way out from this situation is the second 
(coupling) laser using that leads to the electromagneti- 
cally induced transparency (EIT, see more in that case 
in Refs. (HE]) 

Now, let us say a few words about the quantity 
that had been neglected in deriving the equation ([59]) . 
If to take k = (A^i/c), one can find for cesium 
Ek ~ 3, 5 x 10~ 30 eV. So, even at the point Su = 
it is small in comparison with linewidth 7, thus, one 
can note that used approximation is correct. 



6 Conclusion 

Thus, by using the microscopic approach, we studied 
the linear response of the system with bound states 
of particles to disturbing effect of an external electro- 
magnetic field. Our approach was based on a novel 
formulation of the second quantization method in the 
presence of bound states of particles [2] . Such approach 
allowed to obtain the expressions for the macroscopic 
characteristics of the ideal hydrogen-like plasma at low 
temperatures taking into account not only the contri- 
bution of free charged fermions but also their bound 
states (alkali atoms). The expression for the dielec- 
tric permittivity of the ideal gas of alkali atoms in the 
presence of Bose-Einstein condensation phase was also 
obtained. The dispersion equation for the waves prop- 
agating in the system was derived and the existence of 



resonance frequencies was found. 

Our approach gave the opportunity to study the 
propagation properties of the microwave signal, tuned 
up to the transition between two hyperfine ground 
state levels of alkali atoms that are considered in BEC 
state. In contrast to the dipole allowed transitions, it 
was shown that the signal could propagate in such sys- 
tem with rather small energy loss. Such fact allowed 
to introduce the group velocity concept. The slowing 
down conditions for the signal that propagates in BEC 
were studied. Moreover, we revealed the dependence 
of the group velocity sign on the population difference 
sign. 
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